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Abstract 
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^ ' We present a study of bulk thermodynamical quantities in the deconfined phase 

of pure lattice SU{N) gauge theories. We find that the deficit in pressure and 
entropy with respect to their free-gas values, for TV = 4, 8, is remarkably close 
to that of SU(3). This suggests that understanding the strongly interacting 
nature of the deconfined phase, which is crucial for RHIC physics, can be done 
at large N. There, different analytical approaches simplify or become soluble, 
and one can check their predictions and point to their important ingredients. 



> 
(N 



1 Introduction 



Recent discoveries at RHIC point to the inadequacy of using a weakly coupled plasma description for the 
deconfined phase of QCD (see for example Q and references therein). Instead, properties of collective 
flow, can be described very well with ideal hydrodynamics, which suggests that the phase is in fact 
strongly interacting. Related to these are theoretical studies in the framework of lattice QCD, which 
found signs of these strong interactions. These include a large deficit in thermodynamical bulk properties, 
such as pressure, entropy, and energy density, with respect to their free gas values, up to temperatures 
of a few times the transition temperature T c Q l . These deficits were found in the pure SU(3) gauge 
theory, as well as with different types of fermions. The deficit in all these cases is quite similar, and is 
hard to reproduce with standard perturbation theory. This is in contrast to the fermionic addition to the 
pressure at a small, nonzero, chemical potential, which is explained very well by perturbation theory (for 
|^ | example [4| and references thererin). 

The analytical approaches used to study this phenomenon include quasi-particle models, Polyakov 
loop models, variants of perturbation theory, models of loosely bound states, gravity dual models, and 
also perturbative small volume calculations of the large- N gauge theory 2 . Some of these approaches 
simplify or even become soluble at the 't Hooft limit of large- N, which is a free theory in the confined 
phase (at least in the mesonic sector), and consists of only planar diagrams in a perturbative expansion of 
the deconfined phase. An additional important point here is that, provided the chiral limit is taken after 
the large- N limit, the planar pure gauge theory is effectively quenched (for example see 0). This makes 
any calculation in the pure gauge case relevant for the physical case, that includes fermions. 

The simplifications of the planar limit, together with the possibility that the pressure deficit is 
'hidden' in the gluonic sector, and the need to check and constraint the analytical explanations, lead us 
to study bulk thermodynamical properties of SU(4), and SU (8) pure lattice gauge theories and compare 
to the study by Boyd et al. of SU(3) 0. Recent calculations of various properties of SU(N) gauge 
theories 1 6 1— 1 10 1 have demonstrated that N = 8 is very close to N = oo for most purposes. These 
also provide information on the transition coupling for various lattice sizes and N. Thus our calculations 
should provide us with an accurate picture of what happens to bulk thermodynamics at N = oo. For a 
more detailed version of this work, which was conducted with Dr. M. Teper, we refer to [11]. 

'An additional evidence for a 'strongly coupled plasma' is the survival of hadronic states above T c [3]. 
2 An list of some references for these analytical efforts appears in 1111 . 
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2 Lattice setup and methodology 

We define the gauge theory on a discretised periodic Euclidean four dimensional space-time with L 3 x L t 
sites, and perform Monte-Carlo simulations of a simple Wilson action. We use the Kennedy-Pendelton 
heat bath algorithm for the link updates, followed by five over-relaxations of all the SU (2) subgroups of 
SU(N). To evaluate bulk thermodynamics there are several approaches one can use. The "differential 
method" was used in the past, and is known to produce the unphysical result of a negative pressure close 
to T c . The reason for that is the use of a perturbative beta function in a regime where the coupling is not 
small. In our case, larger values of N drive us to smaller values of L t = 5, and therefore to a lattice 
spacing of a ~ 1/(5T C )- too coarse for the differential method. An additional possibility is the use 
of a direct evaluation of the density of states. Unfortunately, modern methods like the Wang-Landau 
algorithm did not converge for the case of SU(8). This leaves us with the "integral method" (see for 
example 1 2 1) in which the pressure p and interaction measure A are 

p/T 4 = 6L 4 fdp'{(u p ) T ~ (u p ) ), A/T 4 = 6Lt((u p (P)} - {u p {f3)) T ) x ■ (1) 

Here (u p )t is the plaquette average on a T > lattice with L t < L s , while (u p )q is measured on a 
lattice with relatively large L t = L s . Let us mention that in the expression for the pressure in Eq. Q we 
dropped the pressure at temperature To corresponding to /3q. To justify this one usually takes (3q much 
lower than the transition coupling j3 c , which means that the systematic error will be very small. 

We performed calculations of (u p }t in SU(4) on 16 3 5 lattices and in SU(8) on 8 3 5 lattices for 
a range of /3 values corresponding to T/T c G [0.89, 1.98] for SU{4), and to T/T c G [0.97, 1.57] for 
SU(8). Since we use L t = 5, while the data for SU(3) in [2] is for L t = 4, 6, 8, we also performed sim- 
ulations for SU(3) on 20 3 5 lattices with T/T c G [1, 2]. The measured plaquette averages are presented 
in fnl . 

We performed the 'T = 0' calculations of (u p )o on 20 4 lattices for SU(3), and on 16 4 lattices 
for SU(4), which are known to be effectively at T = (9]© | for the couplings involved. For SU(8) 
however, using 8 4 lattices would not be adequate for the largest /3-values, and we take instead the SU(8) 
calculations on larger lattices in [ 10], and interpolate between the values of (3 used there with the ansatz 
{u p )o(P) = (u p )i)- T -(P) + ij&(aV°) 4 + C45 8 + c 5 g w , where (u p )fr T -(P) is the lattice perturbative 
result to 0(g 6 ) from (T5| and N = 8. Our best fit has x 2 / dof = °- 93 with dof = 2 > and the best fit 
parameters are C4 = —6.92, C5 = 26.15, and a gluon condensate of = 0.72. 

3 Finite volume effects 

For N = 4, 8, one is able to use lattice volumes much smaller than what one needs for SU(3) (like 
those in \2\) as the longest correlation length decreases rapidly with N QUI. This is also theoretically 
expected, much more generally, as N — > 00 (see for example [5 1). The main remaining concern has to 
do with tunneling configurations, which occur only at T c when V — > 00. On our finite volumes, this is 
no longer true, and we minimise finite-F corrections by calculating the average plaquettes only in field 
configurations that are confining, for T < T c , or deconfining, for T > T c . For SU(3), where the phase 
transition is only weakly first order, it is not practical to attempt to separate phases. This will smear the 
apparent variation of the pressure across T c in the case of SU(3). 

To confirm that our finite volume effects are under control we have compared the SU(8) value of 
(u p ((3)} as measured in the deconfined phase of the our 8 3 x 5 lattice with other L 3 x 5 results from other 
studies llT3l . and find that the results are consistent at the 2 sigma level. We perform a similar check for 
the confined phase on the same lattices, and again find that finite volume effects are small, mostly on 
a one sigma level. A similar check for the confined phase on L 4 lattices leads to conclude that a size 
L = 8 in SU(8) is not large enough, as we find that the plaquette average has a significant change (on 
a 16 sigma level) from L = 8 to L = 16 for our largest value of 0. By contrast, for SU(4) the finite 
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volume effects seem not to be large on the 16 4 lattice as we checked for our largest value of /3 = 11.30. 
There the value of the plaquette on a 20 4 lattice is consistent within ~ 2.3 sigma with the value on a 
16 4 lattice. The data supporting these checks, together with the checks related to the next paragraph, is 
presented in ITU . 

4 Finite lattice spacing corrections 

To evaluate the scaling of the temperature T with the coupling j3 = 2N/g 2 , and the derivative in Eq. Q, 
we use calculations of the string tension, a, in lattice units (e.g. [8|). A more natural scale here is 
the transition temperature T c . However this involves a calculation of the transition coupling j3 c for many 
values of L t , and for several values of L s , which is a very large scale project. Also the difference between 
this scale fixing and the one according to a comes from the 0(a 2 ) corrections of T c / \fa. A check of the 
latter for SU (3), SU (8) shows that T c j \fa for L t = 5, 8, are the same within the errors. Indeed when 
we compare T/T c {j3) from 0, and as determined here for N = 3, we find that they fall on top of each 
other. This is not the case for 5*7(4), but there the 5a difference in T c /y/a of L t = 5, and L t = 8 is 
only at the level of 2%. As a result we may slightly overestimate T/T c when it is ~ 8/5 for this N = 4. 
Finally the value of T c j ' xja for the different gauge groups and L t = 5 in taken from (7]|E). 

In presenting our results for the pressure, we shall normalize to the lattice Stephan-Boltzmann 
result given by (p/T 4 ) free _ gas = (TV 2 — l)|g x Ri(L t ). Here R[ includes the effects of discretization 
errors in the integral method [14 1 1 1. The same normalisation is applied for the interaction measure A, 
and for the energy density e = A + 3p, and entropy density s = (A + 4p)/T. 

To check that we understand the systematics of finite lattice spacing corrections when presenting 
the normalized results, we compare A/T 4 for 577(3) with L t = 4, 6 from 0, L t = 5 from the current 
study (in both these cases T was increased by increasing 0), and L t = 2, 3, 4 from [8] (where T was 
increased by decreasing L t ). The result of the comparison is presented in the left panel of Fig.^ where 
one indeed sees that the normalization by the Stephan-Boltzmann free gas shows a systematic increase 
of A towards the continuum limit. 



■ N=3, [3=5.800 
. N=4,[3=10.635 
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Fig. 1: A for N = 3 (on the left), and TV = 3,4, 8 (on the right) normalized to the lattice Stephan-Boltzmann pressure, 
including the full discretization errors given by Ri. The lines are the 5(7(3) results and Lt = 4, 6 from |2|. 
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N=3, 32 3 x6, Boyd et al. 
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- ■ N=3, 16 3 x4, Boyd et al. 
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5 Results 

Before we present our results for the bulk thermodynamical quantities, let us note that as the ideas of 
reduction at large- N (for example see the recent systematic studies in [5]) have an interesting realization 
here. By relabeling the axes, the authors of IT5ll have argued that at N = oo Wilson loops should not 
change as a function of T for T < T c . This will be reflected in our case in a {u p ) T=T - which is very 
similar to {u p )q. Indeed when we compare the two we find that the difference is on a ~ 2a level for 
577(4), 5*7(8), while from the data of CH it is on the 15cr level for 5*7(3). This also means that the 
systematic error, ignored in Eq. Q is very small for the larger values of N. 
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All the results we present are as a function of T/T c . It is therefore important to note that the 
dependence of T c on N is very weak, indeed like many other quantities in the pure gauge theory 1161 . 
We present our N = 4 and N = 8 results for A/T 4 in the right panel of Fig.[0 and on Fig.|2] In the 
former we change L t so to increase T up to T = 2.5T C , while in the latter T is increased by changing j3. 
In both cases we find modest/small changes between the different groups. Nonetheless in the vicinity of 
T c , the differences are large, presumably because the transition in SU(3) is only weakly first order. 
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Fig. 2: A/T" 1 normalized to the lattice Stephan-Boltzmann pressure, including the full discretization errors given by Ri, for 
N = 3, 4, 8. The line are the SU{3) results and L t = 6 from |2) 



The plot of p/T 4 is presented in the left panel of Fig. |3] We also show our calculations of the 
SU(3) pressure for L t = 5, as well as the L t = 6 calculations from \2\. In the right panel of Fig. [3] 
we present results for the normalized energy density e, and normalized entropy density s. The lines are 
the SU (3) result of [2] with L t = 6. Again we see very little dependence on the gauge group, implying 
very similar curves for N = oo. One can clearly infer that the pressure in the SU (4) and SU (8) cases 
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Fig. 3: The pressure (on the left), and energy and entropy densities (on the right) normalized to the lattice Stephan-Boltzmann 
pressure, including the full discretization errors given by Ri. In the pressure plot, the symbols' vertical sizes are representing 
the largest error bars (which are received for the highest temperature). The lines are the SU(3) results and L t = 6 from [2]. 



is remarkably close to that in SU(3) and hence that the well4tnown pressure deficit observed in SU(3) 
is in fact a property of the large- N planar theory. This implies that the dynamics that drives the decon- 
fined system far from its noninteracting gluon gas limit, must remain equally important in the N = oo 
planar theory. This is encouraging since this limit is simpler to approach analytically, for example using 
gravity duals, and also because it can serve to constraint and point to important ingredients of analyti- 
cal approaches. For example, in perturbation theory, it tells us that the important contributions must be 



5 



planar, and although the current calculated contributions are indeed all planar, this is not guaranteed at 
0{g 7 ) 03 • I n models focusing on resonances and bound states, it must be that the dominant states are 
coloured, since the contribution of colour singlets will vanish as N —> oo. Models using 'quasi-particles' 
should place these in colour representations that do not exclude their presence at N = oo, and in fact 
give them T-dependent properties which depend weakly on N. Also, topological fluctuations should 
play no role in this deficit since the evidence is that there are no topological fluctuations of any size in 
the deconfined phase at large-iV llT8llT9l . 
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